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Abstract. Within the context of loop quantum gravity there are several op- 
erators which measure geometry quantities. This work examines two of these 
operators, volume and angle, to study quantum geometry at a single spin net- 
work vertex - "an atom of geometry." Several aspects of the angle operator 
^ |. are examined in detail including minimum angles, level spacing, and the dis- 

^ ■ tribution of angles. The high spin limit of the volume operator is also studied 

' for monochromatic vertices. The results show that demands of the correct 

scaling relations between area and volume and requirements of the expected 
behavior of angles in three dimensional flat space require high- valence vertices 
with total spins of approximately 10^". 
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lO . 1- Introduction 

O 

. Several promising theories of quantum gravity predict that space is discrete. 

I Within the context of loop quantum gravity, operators for area ||, |j, volume 

|l| ^, |j, and possibly even angle |^ have discrete spectra. Given this radical 
departure from the paradigm of continuous space that has served us so well, it 
O . is important to study consequences of the suggested underlying discrete nature of 

space. This paper reports on a study on the angle and volume operators. More 



5-H . extensive discussion may be found in Ref. Q. Some earlier results were reported 

99' in Ref. §. 

The eigenbasis for all geometric operators is the spin network basis. In this 
k> ■ context spin networks are not merely a graphical notation for representation theory 

. but also contain information that determines the properties of space. All geometric 

' and matter degrees of freedom are contained in the network. Indeed, as traditional 

spin networks are embedded in a three dimensional space knots in the network 
contain physical information about the state of quantum geometry. 

Geometric quantities are concentrated in the edges and vertices of the spin net- 
work. If an edge intersects a surface, the edge endows that surface with area. If a 
vertex (with valence larger than 3) is contained within a region, then the spin net- 
work vertex contributes volume to that region. The angle operator also is defined 
at a single spin network vertex. Since volume and angle are located at vertices, 
the spin network vertex may be called an "atom of geometry." In addition if we 
look beyond the kinematic level, most of the current proposals to implement the 
Hamiltonian constraint contain operator expressions which act at the spin network 
vertex. 

This work is a first step in what we call "single vertex studies." Motivated by 
the fact that volume and angle operators act at a single vertex and the central 
importance of the spin network vertex in dynamics, we studied the relation on 
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Figure 1. Two useful intertwiners. (a) Internal decomposition 
of a vertex operated on by the angle operator; the dotted line 
denotes the boundary between external and internal edges. All of 
the edges Xi,X2, ■ ■ ■ which intersect the surface 5*1 are "collected" 
into the internal edge labelled with ni; similarly, all the edges 
which intersect ^2 are collected into the edge labelled n2, and the 
edges which intersect Sr are collected into the edge labelled rir- In 
the text, we denote si = S2 = Vi^ ^^'^ Sr = X ^i- This is 

called the snowflake basis, (b) The external edges Po, Pi , . . . , Pn-i 
and the internal edges 12,13, ■■ . , in-2 of the comb basis. 

the geometric quantities defined on this space. In such a simple context we can 
answer such questions as: "What are the properties of spin network states which 
approximate continuous space?" and, "What are some of the surprises that arise 
out of this new quantum geometry?" 

In the next section we give the necessary notation, a very brief review of the 
angle operator, and the main analysis. In Section |^ we do the same for the volume 
operator. Section ^ contains a summary of the key points of the study and some 
additional consequences. 

2. Angle Results 

Before reporting on the single vertex study of the angle and volume operators 
we need to introduce some notation on the spin network vertex. The number of 
incident edges is called the valence. While the internal structure of a trivalent vertex 
is unique, higher valence vertices can be constructed in more than one way. These 
different internal structures are described by the intertwiner, which labels the ways 
in which the incident representations are connected. Intertwiners therefore label a 
basis for the vertex. 

The intertwiner may be described by a set of trivalent vertices connected by 
"internal edges." When the edges are to be partitioned into three categories, as 
in the case of the angle operator, we can define a convenient basis in which the 
external edges of one category are collected in a branch which ends in one principal, 
internal edge. The three principal edges are then connected in a vertex we call the 
"intertwiner core." This (class of) intertwiners is called the "snowflake" basis shown 
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Figure 2. A vertex with three partitions of edges identified by 
the surfaces 5*1, 5*2, and Sr- The core of the intertwiner for the 
vertex is chosen so that all edges passing through Si connect to a 
single internal edge (with i = 1, 2, r) as shown in Figure 



in Figure la. The core is the only part of the intertwiner which must be specified 
before completing the diagrammatic calculation of the angle spectrum. When all 
incident edges have the same label the vertex is called "monochromatic." 

2.1. Angle operator: Definition. The idea that one could measure angle from 
the discrete structure of spin networks was first suggested by Penrose This 
work culminated in the Spin Geometry Theorem, which stated that angles in 3- 
dimensional space could be approximated to arbitrary accuracy if the spin network 
was sufiiciently correlated and if the spins were sufficiently large . This result 
may be seen as the "high spin limit" of the observation that, for an EPR pair, we 
know the relative orientation of the particles' spins but nothing about the absolute 
orientation of the two spins in a background space. 

The spin geometry construction carries over into the context of quantum geome- 
try: one can measure the angle between two internal edges of a spin network state. 
At the outset we note that there is more than one definition of the angle operator 
|6|. In this paper we use the quantization that most closely matches the operator 
in the spin geometry theorem. This angle operator is a quantization of the classical 
expression for an angle 6y measured at a point v. The operator assigns an angle to 
two bundles of edges incident to the vertex v. These edges are identified by the two 
surfaces 5*1 and ^2 shown in Figure |. Likewise, there are two angular momenta 
operators Ji and J2 associated to these surfaces. Acting on the spin network state 
\s) the angle is 1^ 

:=arccos (1) 
IJ1IIJ2I 

All edges incident to v are partitioned into three categories corresponding to 
which surface they pass through. In each category, the edges connect to a single 
internal edge in the intertwiner core. These edges are denoted ni, 712, and n^- For 
reasons that will be clear shortly, the spin is known as the "geometric support." 
Using the usual relations for angular momentum operators and an intertwiner core 
labelled by ni, 112 and n^, one immediately finds the spectrum 



?l"^^''|s) = arccos 



Tiririr + 2) — ni{ni + 2) — ^2(712 + 2) 



2^ni(ni + 2) n2(n2 + 2) 



\s). (2) 



We use the integer labels for which ji = ni/2. The key idea of the angle operator 
is to measure the relative spins of internal edges identified with the surfaces Si, S2, 
and Sr- 
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Figure 3. Angle operator spectrum for increasing total vertex 
spin n. The complete spectrum is plotted for total spins from 3 to 
100. 

Note that the quantities rii, n2, and rir are the internal edges that "collect" the 
spins from each of the three surface patches. They are related to the total edge flux 
intersecting each patch by the relations 

ni < si, 712 < S2, and rir < Sr (3) 

where si, S2, and Sr are the respective surface fluxes. 

Normally the discreteness in geometric observables such as area is on the Planck 
scale ip ^ 10~^^ m. However, angle operators have no dependence on the scale of 
the theory £p. The result is purely combinatorial! This result provides a new test 
bed to investigate quantum geometry^] 

Two aspects of the spectrum are immediately obvious from the spectra in Figure 
|[ First, the smallest angle can be quite large; it is "hard" to model small angles. 
Second, the gap between neighboring angles is non-uniform. It is clear from a glance 
that in order to model continuous space both of these effects must be minimized. 
In addition to these more obvious aspects, there is a requirement of these angles as 
embedded in three dimensional space on the distribution of angles. Each of these 
aspects are discussed in the following subsections. 



^ There is an important caveat: Retaining the difTeomorphism invariancc of the classical theory 
in the construction of the state space, there is, geii ericallv. a set of continuous parameters or moduli 
space associated with higher valence vertices . These parameters contain information on the 
embedding of the spin network graph in space. On such a space the angle spectrum is highly 
degenerate. It is by no means clear whether this embedding information is physically relevant 
|12| . Indeed, the state space of quantum gravity may be described by abstract or non-embedded 
graphs. 
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2.2. Small Angles. To obtain a small angle, we wish for the argument of the 
inverse cosine in (H), 

^ nrjur + 2) - ni(ni + 2) - n2{n2 + 2) 

2^ni(ni +2)n2(n2 + 2) 

to be as close to 1 as possible. However, the internal edge spins are subject to the 
triangle inequalities 

2ni, 2^2, 2nr < ni + n2 + fir- 
To maximize Q, then, we wish to maximize its numerator and minimize its de- 
nominator while simultaneously respecting the triangle inequalities. This occurs 
when 

Ur = — and rii = 71,2 « — (5) 

where the "spin sum" n is defined as n — ni + n2 + At these values of ni, n2, 
and Tir, the argument of the inverse cosine then reduces to 

Q = ^. (6) 

n + 8 

For any angle e, then, we have the relation 

e > arccos 



+ 8, 

For small e, cos a; is a uniformly decreasing function; hence, for e <C 1, we have 

4 

e>^_. (7) 

From this relation we see that the minimum possible angle between two patches 
is proportional to n^^/^ for large n. Moreover, we note that if > si + S2, then 
n < 2(si + S2) which implies that for a "macroscopic" vertex, where Sr 2> si,S2, 
the minimum observable angle depends only on the edge flux through the patches 
whose angular separation we are measuring. Since the total spin is bounded from 
above by twice the spin fluxes through the cones, the minimum observable angle 
depends on the area flux through the regions we are measuring. For instance, to 
achieve an angle as small as 10^^" radians, the required flux is roughly 10^*^. 

2.3. Mean Angular Resolution. The second property of the spectrum is the 
non-uniform level spacing. It is possible to quantify the angular resolution of a 
given vertex by counting the number of possible results of an angle measurement. 
We define the mean angular resolution S oi a vertex to be the average separation 
between each possible angle and the next greatest possible angle. In other words, 
3 is the average width of a "gap" between all possible angles in the interval [0,7r]. 
This quantity is not the angular resolution for a given angle measurement; we do 
not compute the angle gap given an angle eigenvalue. 

The simplest way of finding the mean angular resolution for a given vertex is to 
note that an n-valent vertex has n — 2 trivalent vertices in its intertwiner. Since 
each one of these internal vertices can be looked at as a core vertex, we can conclude 
that such a vertex has at most 3(n — 2) possible eigenvalues for the angle operator. 



■^Note that while the value for rir is exact, the values of ni and n2 are only approximate; this 
is because n must be even, but may not necessarily be divisible by 4, and because ni, 712, and Ur 
must be integers. However, for large values of n, the values of ni and 712 that maximize Q will 
approach n/4. 
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Since the 3n— 6 eigenvalues associated with the vertex divide the interval into 3n— 5 
gaps, we conclude that the mean angular separation between the eigenvalues of the 
angle operator for an n-valent vertex is 

3n — 5 

The inequality accounts for possible degeneracies due to the fact that internal ver- 
tices may yield the same angles. 

This method is limited. It only examines the separation between the eigenvalues 
of the angle operators that have a given vertex as their eigenstate. However, other 
angle measurements are possible for a given vertex. For example, if we measure the 
angle between two edges that do not immediately intersect in the vertex's inter- 
twiner, this measurement will not always yield a definite eigenvalue. Instead, there 
will be an expectation value associated with this measurement. The question then 
becomes: for a given vertex, can we find the total number of possible expectation 
values? 

To measure an angle we must partition the edges into three groups (one each for 
Si,S2, and Sr)- In the snowflake basis, we will generally obtain a superposition of 
angle eigenstates. Nonetheless, given a partition we will measure a single, definite 
expectation value, despite the superposition of eigenvalues. We conclude that the 
number of possible expectation values is equal to the number of distinct partitions 
of the edges. 

Since the possible range of colours for the collecting edge is not dependent on 
the order of the edges along the branch, the number of distinct partitions is given 

by 

OO 

^^ = l[Biq,), (8) 

where qi is the number of edges with colour i, and B(n) is the number of distinct 
ways of partitioning n into three distinct bins. Elementary combinatorics tells us 
that the number of ways to do this is 

5(n)=(^"+^)-^)=^!^±l^. (9) 

We must divide the overall product by two, since exchange of the sets corresponding 

to the Til branch and the n2 branch yields the same angle. Thus, the maximum 
number of possible expectation values for a given vertex is 



oo 



+ (10) 

In the monochromatic case, where all the edges have the same colour, Qi is equal 
to the vertex valence n for one value of i (the edge colour) and zero otherwise. The 
formula reduces to 

^^(n+l)(n + 2) 
4 

and the mean angular resolution (between expectation values) is given by 

47r , , 

+ 3n + 3 

where we have again used an inequality because of possible angle degeneracies. 
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It is important not to confuse the "mean angular resolution" as we have defined 
it with the average distance from a random point in the interval [0, n] to the nearest 
angle. What we have found allows us to get an idea of how tightly, on average, the 
eigenvalues of the angle operator for a given vertex are spaced. To find this other 
value (which, if found, would help us to estimate the level to which an arbitrary 
angle can be approximated) we would need 

, , f prob. of being \ /mean separation from\ 

^ \located in the gap / \ either end of gap J ' 



all gaps 

This is equivalent to the expression 



(A) = E S (14) 



a 

where Wa is the width of the gap a. This expression depends heavily on the precise 
location of the eigenvalues in the interval [0, tt], and is therefore difficult to compute 
analytically. 

2.4. Angle Distribution. In the classical continuum model of 3-dimensional space, 
the distribution of solid angles is proportional to sin 9 

7'(6l)d6' = sin6id6i. (15) 

If our quantized angle operator is to be viable, it must reproduce this distribution 
in some "classical limit." To find out whether this is so, we must examine not only 
the location of possible eigenvalues in the interval [0,7r], as we have in the past 
two subsections, but also the likelihood with which these angles occur. Due to the 
complexity of this problem, we will consider only the simplest possible case which 
could still conceivably have a classical limit: that of an n-valent monochromatic 
vertex with edge colour 1. 

We note that every vertex can be transformed into the snowflake basis by re- 
peated application of the recoupling theorem. For a random vertex, however, we 
do not know the exact internal structure a priori; there is no reason to assume that 
any given intertwiner is preferred. We will therefore assume that all intertwincrs 
are equally likely. With this simplification in mind, we can conclude that the prob- 
ability of measuring the angle associated with an intertwiner core ni, n2, is given 

by 

/number of intertwiners\ , 
7'0ni,n2,n, cx . . 16 

\ with core Hi, 77.2, / 

We wish to know, then, how many intertwincrs exist with a given rti, n2, and 
rir- It can be shown that if i is the number of edges (of colour 1) entering a branch, 
the number of distinct branches {02, as, . . . , a^} that end in a given ni — a; is Q 

«M.)^^(4,^;,) (1.) 

To turn this into a normalized probability distribution, we approximate the bino- 
mial coefficient in the above expression by a Gaussian distribution multiplied by a 
normalization factor 
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where A is a normalization factor. Applying this to (|T^, we have 



fli + 1 



2(i + 1) 



A has been subsumed into the exponential term. 

Using this probability, we can numerically compute the distribution of the an- 
gles in the interval [0, tt]. However, the sheer number of possible vertices for a 
macroscopic vertex prohibits calculating every angle associated with such a vertex. 
Instead, we sampled the possible angles randomly, resulting in a certain amount 
of scatter in the data points. This allowed us to gain a qualitative understanding 
of the distribution. The results for various values of Si, S2, and Sr are shown in 
Figures and (||). 

Figure (|^) shows the distribution for si = S2 — Sr- This distribution certainly 
does not match the expected classical distribution; most notable is the fact that 
the distribution peaks at approximately 2.5 rad « 140°. Figure (||) shows the dis- 
tribution for si < S2 <SC Sr- This distribution has surprisingly good correspondence 
with the expected classical angle distribution. 

These results, while not conclusive, are certainly indicative that the angle opera- 
tor can indeed reproduce the classical angle distribution in the case where Si, S2 ^ 1 
and Si,S2 <C Sr- This first condition corresponds to what normally thinks of as a 
"classical limit," well away from the regions where quantum effects dominate. The 
second condition can be thought of as a requirement for a sufficient amount of 
"background geometry": not only do we need to make sure that the angles we 
are measuring include a large amount of spin, but we also need to ensure that the 
background spacetime upon which we measure the angles is sufficiently classical. 



3. Volume Results 

3.1. Volume Operator: Definition. After briefly providing the definition of the 
Rovelli-Smolin-DePietri volume operator, we present the large spin results. For 
more on the operator itself, see Ref. . 

The volume operator acts on a vertex. In this case, however, it is more convenient 
to use the basis pictured in Figure |l|b, where all the edges are lined up sequentially. 
One also defines the volume-squared operator (known as W) first. The volume 
eigenvalues are the square roots of the absolute values of the eigenvalues of W . 

For a vertex with n edges labelled {Po,Pi,... ,P„_i}, we define an operator 
W^j-st] for each triple of edges such that 0<r<s<i<n — 1. W\^rst] grasps edges 
r, s, and t: 



Pq p, p, p. 



Po p. P, P, 



(19) 



'1 '4 



2 2 
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Figure 4. Angle distribution for si ~ S2 = s^. 
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We can write the left hand side of this equation as a superposition of other vertices 
with the same internal decomposition as those in Figure |^b 

p« p, p, p />„_, p« p, p p. p., 



k2,... ,k„ 



A', 



(20) 

This is essentially a matrix equation with the entries of the matrix indexed by the 
internal intertwiner edges. These entries turn out to be Q 

k2 Pt h 



w, 



(n) k2---kn-2 

[rst]i2---'f-n-2 



-PrPsPti ^2 Pt *3 yA^f<5f;...<": 



r Pr 


Pr 






Ps k4 \ 


\ k2 


«2 






«3 2 / 



(21) 



eik2,k3,Pt) 

In this formula we have used the 9-j symbol, which is equivalent to the spin network 



k2 Pt k3 
«2 Pt is 
2 2 2 



Liz I r h I j 

_\ i 



(22) 



The eigenvalues of the volume-squared operator are then proportional to the 
sum of the absolute values of the eigenvalues of the W^-operators 



E 



0<r<s<t<n-l 



16 



rst] 



(23) 



and the volume eigenvalues are the square roots of the eigenvalues of the V oper- 
ator. 

3.2. The T^/'-matrix and Eigenvalue Bounds. One can see that the expression 
for the volume of Eq. (|2^) is horrendously difhcult to handle in all its generality, 
for three main reasons. First, for an n-valent vertex with arbitrary edge colours, 
there are (g) possible choices for r, s, and t; each of these choices could conceivably 
produce a different M^/"!] . Second, the eigenvalues of the V'^ operator are difficult 



to find; even those of the M^-matrix are distinctly non-trivial, since all that Eq. (21) 
gives us is the matrix entries. Third, the matrix entries themselves, as given by Eq. 
(pT|), are not easy to calculate analytically. 

With these problems, it is to our advantage to make some simplifications. Instead 
of considering an arbitrary vertex, we will examine the case of a monochromatic 
vertex for which all of the matrices W^j^r"^] are identical (since our "graspings" are 
always on edges of the same colour) The formula for volume reduces to 



y2 








n(n — l)(n — 2) 




16 


96 



(24) 



•^In this case, we will often suppress the indices [rst], since they are always the same; in other 
words, when referring to the monochromatic case, we will define VF'"' = ■ 

' ' [rst\ 
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This solves our first problem. The second problem, however, is somewhat less 
tractable. The matrix has a row for each possible intertwiner core - usually a 

fairly large number. To get around this, we will examine bounds on the eigenvalues 
instead of the eigenvalues themselves. For any n x n matrix with entries a^ , any 
given eigenvalue A of this matrix satisfies 

n n 

|A|<max> \aij\ and |A| < max > \aij\- (25) 
]=i 1=1 

In other words, the magnitude of every eigenvalue must be less than the sum of the 
absolute values of some column and of some row. Using this, we can find an upper 
bound M on the magnitudes of the eigenvalues of W'^"'\ In the monochromatic 
case the eigenvalues Xwa of the iy(") operator are related to the eigenvalues Ay2„ 
of the V"^ operator by 

n(n-l)(n-2) ,^ 

^V^a~ \^Wa\ i^Oj 

as can be seen from (|4|). Hence, we can place a limit on the eigenvalues of the 
volume operator, Xva ~ V \v^a. 



/n(n- l)(n-2) 
A^.<^^ ^ '-M. (27) 

The determination of the upper bound M is also rather complex. Simpler expres- 
sions do exist for the case of a 4-valent vertex; in this case, we can easily calculate 
the matrix elements explicitly. For higher valence vertices, numerical techniques 
come in handy. 

3.3. Volume Eigenvalue Bounds for 4-valent Vertices. In the case of a 4- 
valent vertex, De Pietri has shown that the entries of the W'^^'^ matrix for a 
4-valent vertex with edge colours a, h, c, and d are given explicitly by the formula 

^(4) _ -,(-l)(a+'>+c+^)/ 2 



[•'^'l*-^ 32^^(1+2) 

[[a + b + t + 'i)[c + d + t + 'i)[l + a + b-t) 
{l + c + d-t){l + a + t-b)[l + b + t-a) 
{l + c + t-d){l + d + t~c)Y . (28) 

where t + e and t — e correspond to the internal edge of the vertex. It can be shown 
that these elements are zero unless e = ±1; thus, each row or column will have at 
most two non-zero elements. To find M in this case all we have to do is to find the 
maximum value of this polynomial. 

In the monochromatic case, a = 6 = c = d; we call this edge colour m. The 
polynomial in question becomes 

(2m + t + 3)(2m - t + l)(t + 1)^ 



TT.(4) t + , 
[012] t-f 



32^i(i + 2) 



(29) 



Thinking of as a continuous function W{t)^ this expression is maximized at 



y (2m + 2)2 + 4 ^ l^(2m + 2)^ - 16(2m + 2)^TT6, (30) 
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Figure 6. The eigenvalue limit obtained in ( |29D and (|30|), com- 
pared with the maximum volume eigenvalue for 1 < m < 100. The 
lower dotted line is the result of a curve fit of the maximum eigen- 
values to a curve of the form A = km?/'^^ where k is the only free 
parameter of the fit. This fit yields a value of fc = 0.22155. 



and the maximum value is, of course, W{tr 
m ^ 1, this maximum value scales as 



\Wrr 



6\/3' 



In the large-spin limit, where 



(31) 



The maximum absolute row or column sum (i.e. the quantity M) is at most twice 
this number; therefore, ( p7| ) tells us that for a 4-valent monochromatic vertex with 
edge colour m, any eigenvalue \v of the volume operator on this vertex will satisfy 
(more or less) the inequality 

\y<ll^=. (32) 

This is a very encouraging result when we consider it in relation to the area 
operator. The eigenvalues of the area operator are given by ||, |3| 



2) 



(33) 



where the values rrii are all the intersections of edges with the surface in question. 
For a 4-valent vertex with large edge colour m, we have 

Xa oc m. (34) 

Together with (^), this implies that the volume scales no faster than A'^/^, which 
is what we would expect from classical spatial geometry. 

We must ask ourselves, however, how "good" this bound is: do the eigenvalues 
actually scale proportionally to m'^/^, and if so, is the proportionality coefficient 
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we have found (i.e. 1/v 12-\/3) close to the actual proportionality coefficient? Fig- 
ure H shows the largest eigenvalue of the volume operator (found by numerical 
calculation) compared to the bound stated in ( p9| ) and (^0[). As we can see, the 
bound fits the volume eigenvalues very tightly, diverging significantly from them for 
71 < fO and becoming practically indistinguishable from the maximum eigenvalue 
for n > 50. 



3.4. Eigenvalue Bounds for 7i-Valent Vertices. We now turn our attention 
to the more general case of the ri-valcnt monochromatic vertex. In this case, we 
cannot use the simple formula of the four-valent case, but must instead use the 
more general form. There are some simplifications that can be made which aid us 
in computing these matrix entries and their bounds. Numerical methods prove to 
be fruitful. 

We recall that we are examining the special case where Pr = Ps = Pt = Po = m. 
We note that the expression for the 9-j symbol in Eq. ( pl| ) is fairly symmetric, and 
can be simplified. This symbol corresponds to the spin network 



r k2 


Pt 
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\ 

P, 




Pt 




v. 
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Si 






'3 U 










\. 




J 



(35) 



Using usual recoupling moves we can move the graspings of the 2-edges and close 
the resulting "triangles" to obtain 
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Plugging this back in, (|2l|) becomes 



where 
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(36) 



We can write the Tets in this equation in terms of Wigner 6-j symbols. At least 
one of the arguments of the resulting Wigner Q-j symbols will be 1, and analytic 
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formulas for this case have been compiled (see, for example, Varshalovich |p^). 
Equation ( pT] ) then becomes 



where Q is given in (|3^), and TZ and S are given by 

r m/2 /C2/2 fc3/2 \ / m/2 m/2 fca/Z \ / ki/2 m/2 /C3/2 , 
1 1 Z3/2 ^2/2 ]^ \ 1 ^2/2 m/2 \ 1 Z3/2 m/2 



c /J/ o\ / S{^2, is, m)e{m,m, 12)0(13, 14, m) 

S = e[m,m,2)J — — — — -. 39 

Note that we are using a W subscript to denote where we have used Wigner 6-j 
symbols instead of the usual spin-network 6-j symbols; see Ref. for the exact 
correspondence between these sets of symbols. 

Finally, we note that there are certain "selection rules" for finding the non-zero 



elements of W in this basis. The presence of a series of Kronecker deltas in (|2l|) 
gives us the obvious selection rules that 

fc4 = 14, fcs^is, •■■ = jn-2- (40) 

It is also evident from the spin diagram in (|2|) that for the 9-j symbol in (^ ) to 
be non-zero, we must have 

fc2-«2 = 0,±2 and A;3 - ^3 = 0, ±2 (41) 

in order to satisfy the triangle inequalities. Finally, it can be shown that 

(_l)a+i.£i±2 b^a~2 



-Tet 



a a b 
2 2 2 




0(a,6,2) 

Hence, if k2 = 12 and fcs = 13, the two terms in parentheses in (|3(1) will be be equal 



b ~ a 
b^a + 2 
otherwise 



(42) 



and therefore cancel out, making the quantity Q zero. We can then rewrite (41) as 

(43) 



A:2-i2 = 0,±2, fc3-«3 = 0,±2, with 
^2 — *2 and ^3 — is are not both zero. 



We conclude that for a given column of the W matrix (i.e. fixed fc2, k^, . . . , kn-2), 
there will be at most eight non-zero elements: 

«2 = ^2 - 2 l2 = fc2 - 2 12 = k2-2 

is. = kz-2 i3 = ks = k^ + 2 



(44) 



42 = k2 l2 = k2 

^3 = fc3 - 2 = fc3 + 2 

«2 = ^2 + 2 12 = k2 + 2 12 = fc2 + 2 

h = k3-2 is = ks is = ks + 2 

The quantities Q, TZ, and S are fairly simple to calculate for a given set of 
arguments. Finding the maximum product of all three, over all possible values 
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of fc2, /ca, A:4, *3, and 14, is a much more daunting taskjj This task is further 
comphcated by the fact that both the quantities Q and TZ have shghtly different 
analytic forms for each of the eight entries in a given column — one form for 
k2 ~ 12 — "2 and — is — 2, one form for fc2 ~ ^2 = 2 and — = 0, and so forth 
(see Varshalovich JTst-) 

Despite this difficulty, we can get an idea of the scaling properties of the absolute 
row and column sums by numerical analysis. By a happy coincidence, we do not 
even need to look at all of the rows and columns of the matrix — only a small subset 
of them. This is because (as mentioned above) the entries of the W^-matrix are only 
dependent on the six values {A;2, ^3, ^4, ^2, 13, «4}; all other pairs of corresponding 
indices must be equal, or the matrix entry will be zero. These facts imply that the 
VF-matrix must be of block-diagonal form, and, more importantly, that each one 
of the blocks is identical. We see that we need only compute the absolute row and 
column sums for the (much smaller) block matrix, instead of the full W^-matrix. 
This is still a non-trivial calculation, as the size of this matrix scales fairly quickly 
(for m = 10 this matrix has already grown to 891 x 891) but it does simplify the 
problem greatly. 

Using this trick, combined with the selection rules in ( ^ ) and (^3|), the bounding 



quantity M (as defined above in section |3.3D may be calculated. Two methods 
were used. The first was the absolute row and column sum discussed previously. 
The second was MATLAB's internal normest function, which estimates the largest 
eigenvalue of a sparse matrix.^ These results are shown in Figure 0. Each of these 
data sets is fitted to a power law relation. We see that for both of these data sets, 
we have an approximate relation M cx m^. 

There is good reason to believe that this power dependence should tend towards 
(exactly) as m — > 00: if it does, then we have the relation 



^^n(n-l)(n-2) 



96 



CVrnFn^ 



= C(total spin)^/^ (45) 

where C is a proportionality constant. As noted in the previous section, this is 
the dependence expected from classical spatial geometry. While our data are not 
conclusive evidence of this relation, they are highly suggestive that this relation 
holds for a monochromatic vertex of any valence. 



4. Conclusions 

We have seen some interesting consequences stemming from the discreteness of 
the angle operator. These results will be referred to as the small angle property, the 
angular resolution property, and the angular distribution property. We now wish 
to assign some scaling behavior to these predicted phenomena; in other words, on 
what length scales do we expect these phenomena to be observable? 

* Note that, in a small mercy granted to us by equation the entries of this form of the 

W-matrix depend only on these six intertwiner strands (and the edge colour m) and not on the 
entire set. 

^ For the mathematically inclined reader, we note that the absolute row and column sums are 
known as the 1-norm and the 00-norm of a given matrix, respectively; that which normest finds 
is the 2-norm. 
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« — "Normest" matrix norm 
^ Abs. row/col. sum 




1 10 

Edge colour 

Figure 7. The eigenvalue limits found by the absolute row and 
column sums and by MATLAB's internal normest function. The 
curve fits are for a power law relation. For the absolute row/column 
sum method, the best fit is M = 0.90432 • x^-^i^e. Jq^. ^^^^ normest 
method, it is M = 0.7f884 • x'^-^^^'^. 



The small angle property is the best-quantified of these three properties. We 
have the equation 

4 4 

e > , > , (46) 

where st is the total surface flux. From our numerical studies of the volume oper- 
ator, we also have the result that 

y cx cx (47) 

If we combine these results then we obtain the result that for large st, the minimum 
observable angle is related very simply to the length scale observed 

e cx cx s^^/^ (48) 

where L is the characteristic length (in Planck units) defined by the volume or 
area. Thus, an accuracy of 10"^" radians corresponds to spin fluxes of 10^" and 
spin network vertices of L ~ 10~^^m. Moreover for large "background geometries" 
the primary limit is not from st but from (si -I- S2) 

ecx (si + S2)"^^^. 

It is also possible to directly relate the minimum angle to the volume contribution 
of the spin network vertex. Using the same relations as above the dependence of 
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minimum angle and volume on the total spin flux implies the "minimum-angle 
volume bound" given by 

> 1- (49) 

One way to read this is that the smaller the minimum angle desired, the larger the 
volume contribution of the vertex, the larger the "quanta of volume." Thus, this 
particular form of the angle operator suggests a very different description of our flat 
space. The effectively continuous spin network state would be a highly ffocculcnt 
network of high valent vertices. 

We note that this "bound" is based only on the considerations of the dependence 
of the spectra on the spin flux. It is not derived by the non-commuting properties 
of the operators. In addition to this bound, we expect that there would be a 
volume-angle unc ertainty relation as there is in the case of the area operators p4| . 



In Section 2.3 we found results for the mean angular resolution associated with a 
vertex. The most important of these were the average spacing between expectation 
values for momochromatic vertices, given by 

Combining these results with the result in Eq. allows us to associate a scaling 
behavior with the mean angular separation and resolution. For large n, the angular 
separation scales as 

(5 cx oc L"*. (51) 

This scaling property of the angular resolution is somewhat striking, since it scales 
proportionally to such a high value of the total spin (i.e. inversely proportional to 
the fourth power of the length.) 

We found that we could reproduce the classical angle distribution of V{9) dO = 
sin 9 d9 under the condition that 

1 < {si, S2} <C Sr- 

In the absence of good analytical control on the precise behavior, it is not possible 
to determine the length scale to this condition. The importance of our angle distri- 
bution results is more clear in light of the small-angle results. The minimum angle 
observable depends on st if si,S2 ~ s^, but on si + S2 if si,S2 <C Sr- To obtain 
the classical distribution the latter case holds and so the sum of the fluxes si and 
S2 must also be large. Thus, to approximate classical angles the total spin and the 
spin through surface patches 5*1 and ^2 must be quite large. 

Finally, in reflecting on these results it is important to keep in mind the nature 
of angle operator under consideration. Classically an angle is defined at the inter- 
section of two lines - a point. In loop quantum gravity it is possible to define the 
angle at a single spin network vertex. As the above analysis shows, these "atoms 
of geometry" are small objects. It seems likely that to model familiar interactions 
which have "angle dependence" it will be necessary to average over many such 
vertices. If this in fact proves to be the case, then the above estimates must be 
modified. 
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